The Effect of Eccentricity on Searches for Gravitational- Waves from Coalescing 
Compact Binaries in Ground-based Detectors 
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Inspiralling compact binaries are expected to circularize before their gravitational-wave signals 
reach the sensitive frequency band of ground-based detectors. Current searches for gravitational 
waves from compact binaries using the LIGO and Virgo detectors therefore use circular templates 
to construct matched filters. Binary formation models have been proposed which suggest that 
some systems detectable by the LIGO-Virgo network may have non-negligible eccentricity. We 
investigate the ability of the restricted 3.5 post-Newtonian order TaylorF2 template bank, used by 
LIGO and Virgo to search for gravitational waves from compact binaries with masses M < 35Mq , to 
detect binaries with non-zero eccentricity. We model the gravitational waves from eccentric binaries 
using the x-model post-Newtonian formalism proposed by Hinder et al. [I. Hinder, F. Hermann, P. 
Laguna, and D. Shoemaker, arXiv:0806.1037vl]. We find that small residual eccentricities (eo < 0.05 
at 40 Hz) do not significantly affect the ability of current LIGO searches to detect gravitational waves 
from coalescing compact binaries with total mass 2Mq < M < 15M0. For eccentricities eo > 0.1, 
the loss in matched filter signal-to-noise ratio due to eccentricity can be significant and so templates 
which include eccentric effects will be required to perform optimal searches for such systems. 

PACS numbers: 04.30.-w, 04.25. Nx. 04.30. Db, 04.80.Nn 
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I. INTRODUCTION 



In October 2007 the Laser Interferometer 
Gravitational- Wave Observatory (LIGO) [l| com- 
pleted its fifth science run, collecting one year of 
coincident data at design sensitivity After a period 
of detector upgrades 0], LIGO and the French- Italian 
Virgo detector Q began a new observing campaign 
in July 2009 and construction of the Advanced LIGO 
(AdvLIGO) detectors [sj is underway. The inspiral and 
coalescence of compact binary systems consisting of 
neutron stars (NS) and/or black holes (BH) is a promis- 
ing source of gravitational waves for the LIGO-Virgo 
network 0]. Optimal searches for inspiral signals 
use the method of matched- filtering [1, @ , in which the 
time-series of the detector is cross-correlated against a 
set of theoretical template waveforms. The utility of 
this technique rests on how well the template waveforms 
model the signal being sought. Significant progress 
has been made in modeling gravitational waves from 
compact binaries using post-Newtonian (PN) theory [l^ 
and numerical relativity (NR) [llj. For binaries that 
evolved through typical main sequence evolution [T^j [13 , 
radiation reaction will cause them to circularize by the 
time the frequency of their gravitational waves enters 
the sensitive frequency band of ground-based detectors 
such as LIGO and Virgo [3, Alternative forma- 

tion mechanisms have been proposed that may yield 
binaries with non-negligible eccentricity in the LIGO 
band [l^ The scattering of stellar mass BHs in 

galactic cores containing a super-massive BH may lead 
to binary formations with high eccentricities; ~ 90% 
with Co > 0.9, where eo denotes the eccentricity of 
the binary when it enters the sensitive band of the 
detector [Ig. In such an encounter, the two BHs radiate 



enough energy to form a bound system, but do not fully 
radiate their eccentricity in GWs until they coalesce. 
Estimates suggest that the rate of such BH binary 
coalescences detectable by AdvLIGO may be as high as 
100 yr~^ (l6j . Binary-binary interactions in the central 
regions of globular clusters may also be a source of 
gravitational waves with non-zero eccentricities. Multi- 
body interactions in globular clusters may result in the 
formation of a stable hierarchical triple. If the orbital 
planes of the inner and outer binary with respect to the 
center of mass are highly inclined with respect to one 
another, Kozai resonance increases the eccentricity of 
the inner binary It has been estimated that ^ 30% 
of binaries formed in this way will have eccentricities 
eo > 0.1 when they enter the AdvLIGO frequency band 
(~ 10 Hz). However, most template families, including 
those used by LIGO to search for inspiral signals, make 
use of the quasi- circular approximation. In this scheme, 
the binary system evolves secularly through a sequence 
of circular orbits 0, [l^; eccentricity is neglected. 

In recent years there has been significant progress in 
modeling eccentric binaries using post-Newtonian the- 
ory. These calculations use a combiiiation of multi-scale 
methods and variation of constants i2Q|, 122, i23| to include 
conservative post-Newtonian effects in the phase evolu- 
tion [l^] in addition to the effect of radiation reaction. 
The conservative dynamics account for small oscillations 
in phase evolution which occur on time-scales of the order 
of an orbital period. At present, the 3 PN conservative 
contributions to the phase evolution have been derived in 
the quasi-Keplerian parameterization [245 ^^'^ the 3 PN 
reactive contribution is expected shortly now that the 
calculation of the 3 PN flux 25] is complete. 

The first investigation of the effect of eccentricity on 
searches for gravitational waves was performed by Mar- 
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tel and Poisson \2^. In their study, they examined the 
effect of the leading-order radiation effects and con- 
cluded that circular templates were effective at capturing 
sources with small residual eccentricity. This problem 
was recently revisited by Cokelaer and Pathak ^271] using 
a discrete template bank in an effort to more accurately 
model an inspiral search. The qualitative conclusions of 
Ref. [2^ were upheld by Ref. [23 , despite small quantita- 
tive differences due to a difference in the numerical pre- 
cision used in the computations [27| . These studies used 
waveforms computed at leading order; both conservative 
effects on the waveform evolution and higher order post- 
Newtonian corrections to the phase evolution were omit- 
ted. Recent work by Tessmer and Gopakumar, based 
on templates which include conservative post-Newtonian 
effects, suggests that circular templates may not be ade- 
quate to detect systems with residual eccentricity j2^. 

In this paper, we investigate whether the waveforms 
used in current LIGO inspiral searches are adequate to 
detect gravitational waves from compact binaries with 
residual eccentricity. To model eccentric signals, we use 
the post-Newtonian eccentric waveform family proposed 
by Hinder et al. [HI, which has been calibrated against 
numerical simulations of eccentric BH binaries. In Sec.HIl 
we review the waveforms used as eccentric signals and 
in Sec. IIIII we review the data analysis techniques used. 
Sec. |IV] presents our results for the current and AdvLIGO 
detectors. 



II. ECCENTRIC POST-NEWTONIAN 
WAVEFORMS 

The effect of gravitational radiation on the orbital mo- 
tion of eccentric binaries was first calculated by Peters 
and Mathews [13, [iBl- Post- Newtonian corrections to 
the Peters and Mathews result were first calculated in 
Ref. and later refined using an improved method of 
variation of constants in Ref. The PN equations of 
motion in Ref. [2^ were written in terms of the mean 
motion n and time eccentricity et and the 3.5 PN phase 
evolution has been calculated using these variables ^3} . 
The recent progress in the numerical evolution of BH 
binaries now allows post-Newtonian models to be cal- 
ibrated against the gravitational waves extracted from 
simulations. Hinder et. al. investigated the agreement 
of eccentric post-Newtonian waveforms with the gravi- 
tational waves extracted from an equal mass BH binary 
with eccentricity eo = 0.1 [2^. Based on this study they 
proposed a modification of the eccentric PN waveforms 
(called the x-model) in which the dynamical quantities 
are written in terms of a variable related to the orbital 
frequency x — (Maj)^/"^ and the time eccentricity et- This 
waveform gives better agreement to the numerical data 
than the waveform written in terms of n and et and so 
we use it to model gravitational waves from eccentric bi- 
naries in our analysis. In the remainder of this section, 
we review the x-model formalism introduced in Ref. 12911 . 



A. Conservative post-Newtonian Dynamics 

The post-Newtonian model of Ref. [2^ is more eas- 
ily understood if we first consider the evolution of two 
point masses mi and m2 in the purely Newtonian case. 
In the Newtonian case, conservation of energy E and an- 
gular momentum J dictate the evolution of the orbital 
elements. We define the mean motion n in terms of the 
orbital period P, n = For Newtonian trajectories 
on an ellipse (0 < eo < 1), n is simply -y/ M/aP' where 
M = mi + 7712 a-nd a is the semi-major axis^. To param- 
eterize the equations we introduce the eccentric anomaly 
u. Written in terms of u, the relative orbital separation 
r and the angular frequency </> take the form 



r = a [1 — ecosu] 
• ny/l — e^ 



[1 



ecosM 



(1) 
(2) 



To complete the dynamical system we relate the eccentric 
anomaly to the mean motion via the Classical Kepler 
Equation 



I = 



esmu, 



(3) 



where the mean anomaly I is given by integrating dl — 
n dt. In the absence of radiation reaction, 77 is time- 
independent and the mean anomaly is simply I = n{t — 
to). These equations form the Keplerian parameteri- 
zation of a Newtonian orbit. Given initial conditions 
(fiQ = (j){to) and Iq = l{to) along with numerical values 
for eo and u we can compute the trajectories of the two 
particles at any time t > to by root finding for u. 

The post-Newtonian dynamics of a system are conser- 
vative if they the respect energy and angular momen- 
tum conservation^. Using the relations given in Ref. [l^l 
for 77 and et in terms of E and J, one can derive 3 PN 
quasi-Keplerian equations of motion by introducing three 
eccentricities, et, e^ and e^, which account for the varia- 
tions in the t, r, and directions. The addition of higher 
order conservative PN corrections lead to periastron pre- 
cession. We let A(j> represent the angle of precession dur- 
ing one period P. The angle swept out over consecutive 
periastron passages is defined 



2% + A(f) 
P 



(4) 



[25l . |29| . |30| . We note that w is a constant in the absence 
of radiation reaction. 

The evolution equations in the conservative dynamics 



^ We work in units with G = c = 1. 

^ This definition fails at 3.5 PN due to center of mass "recoil". A 
more rigorous definition is given in Ref. [221. 
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take the abbreviated form (2£ 



r 

M 



I 

Mi 



— (1 — et cos u)x ^ + riPN + r2Pi^x 

+ rsPNX^ +0{x''), (5) 



(1 — et cosm) 

u — et sinw + hpNX^ + hpNX^ 
Mn = T^/^ -I- n,T^^TT^/'^ 



X 
^9/2 



+ niPNa; 



n2PNa; 



0(x^ 

7/2 



(6) 
(7) 

(8) 



Detailed expressions for the PN coefRcients are given in 
the appendix of Rcf. [29]. We note that the PN co- 
efficients rpN, 0PN, ■ ■ ■ are functions of both et and u, 
whereas n depends on et alone. Conservative PN trajec- 
tories are obtained by first integrating Eq. ([8]) for I, and 
then numerically root solving Eq. ([7]) for u at each time 
step t. 



B. Radiative post- Newtonian Dynamics 

Radiative post-Newtonian dynamics are needed to de- 
scribe evolutions over time-scales in which angular mo- 
mentum and energy are carried away from the system 
by gravitational radiation. The time dependence of E 
and J implies that x and et are no longer integrals of 
motion. The time variation of n and et leads to a secu- 
lar evolution of the orbital parameters. To describe the 
secular evolution of the orbital elements (f), 4>, r, f, and I, 
the equations of motion for x and et must be extended 
into the non-conservative regime. The 2 PN equations 
describing the radiative dynamics read"^ 



Mx 



Met = 



2r] 



jj^ (96 292e2 + 37ef ) x^ ^ xipnx^ 

2 A / 



15 (1 - ei 

il.5PNxl3/2 ^ ^2p^2.7 ^ 0(3.15/2) 



(9) 



(304-Hl21e?) x'^-HeiPNX^ 



15(l-e?f/2 
+ ei.5PNa;"/2 + e2PNa;^ + 0(a;"/2), (10) 

where r\ = /u/M = mim2/{mi + 1712)'^ is the symmetric 
mass ratio. Again, we refer to the appendix of Ref. [2^ 
where the PN coefficients are written out explicitly. In 
the adiabatic approximation, the PN coefficients are in- 
dependent of u, and are solved independently of the Ke- 
pler Equation. We use a fourth order Runge-Kutta- 
Fehlberg (RKF-45) with adaptive step-size control [3l| 
to numerically solve for the radiative dynamics at each 



time t > to. Once x{t) and et{t) are obtained, we sub- 
stitute their values into Eq. ^ and solve the ordinary 
differential equation (ODE) for l{t). The value of l{t) is 
then equated to the RHS of Eq. (O, giving us our PN 
Kepler Equation. We then solve the PN Kepler Equation 
using the Mikkola root finding method [321 ■ The values 
of X, et, and u are then substituted into the conservative 
equations for r and 0. The value of r is obtained numer- 
ically using a five-point stencil method and the phase 
is numerically integrated using the RKF-45. Repeating 
this process, we obtain the dynamics (r, r, (j), 0) at each 
time t > Iq. 



Formulation of the Dynamics in The Zero 
Eccentricity Limit 



In the zero eccentricity limit, the cc-model phasing for- 
malism reduces to a dynamical system closely resembling 
those of the Taylor time-domain approximants. The sys- 
tem of ordinary differential equations reduces to 



Mx 



Me 
M<i) 

Mi 



^-fxHl 



336 4 
34103r/ 13 661ry2 59r;3 



18144 



2016 



18 



= 0, 

= x^/^ _^ 3^5/2 _^ 1 (2877 - 18)a;^/2 



UStt'^t] 48177 27 
32 ^ ~ 2 



^9/2 



(11) 
(12) 

(13) 



(14) 



The evolution equations of the system in the zero eccen- 
tricity limit are very similar to the TaylorT4 approxi- 
mant. The equations which govern the phase evolution, 
(Eqs. (III|),(Il31)) are equal to the 2 PN TaylorT4 equa- 
tions j33j . The eccentric x- model formalism differs in 
that it has an amplitude contribution entering in through 
Eq. We explore the similarity between the two mod- 
els in greater depth in Sec. IIVI where we discuss overlap 
calculations. 



D. Eccentric Binary Waveforms 



^ It is conventional in the gravitational-wave literature to define 
the radiative PN order relative the order at which radiation re- 
action occurs. 



We now discuss the eccentric waveforms used for sig- 
nals in the template bank simulations. The leading order 
time domain gravitational-wave polarizations are given 
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by the following expressions: 



-Mr] 
R 

2rf(j) sin 20 
-2Mr] 



(1 + cos\) 



M 



+ r'^cjp- - 1 cos 20 



M 



2 12 -2 

r — r 



sin^ L > ,(15a) 



cos L<, \ — + r'^cfP' — I sin 20 



-2rr0cos 20 > , 



(15b) 



where i? is the radial distance of the binary, and l is 
the inclination angle of the orbital plane of the binary 
measured from the line of sight [2l|, [sj] . When the orbit 
is circular, final term in Eq. (|15ap vanishes. This can be 
seen from the fact that the radial velocity r is zero and 
the Newtonian equation for the centripetal acceleration 
is M/r = r^0^. For circular motion the polarizations 
simplify to 

— 2M?7r^O^ 

h+ = :k (1 + cos^t)cos20, (16) 



R 

Aix = — COS/, sin 20, 

R 



(17) 



where r is the binary separation and fl is the orbital 
frequency. Notice that we are using the dominant har- 
monic only. Since the post-Newtonian waveforms used 
here do not capture the merger or ringdown of the sig- 
nals, we terminate the waveforms at the Schwarzschild 
innermost stable circular orbit (ISCO) frequency, /isco = 
1/ (6-\/67rM). The waveform h{t) observed at the detec- 
tor is a linear combination of the -|- and x polarizations: 



h{t) = F+i9, if, ^)h+it; L, 00 ) 
Fx(0, V)/ix(i; 00 ) 



(18) 



where F+ and Fx arc the beam-pattern factors of the de- 
tector and 00 is the azimuthal Euler angle of the source. 
For ground based interferometers, the beam factors are 
expressed in terms of the Euler angles of the detector 
[9, ip) and the Euler angle of thepolarization plane ip. 
Following the convention of Ref . [35| , the expressions for 
the beam-pattern factors are 



1 



1 + cos^^) cos 2(f cos 2ip 



cos 6* sin 2(p sin 2-0, 
1 



(p, 0) = - (1 + cos 6) cos 2(f sin 2-0 
— cos 9 cos 2lp cos 2-0. 



(19a) 
(19b) 



In computing the waveforms, we set the angles such that 
the binary is optimally oriented for the -I- polarization 



mode; i.e, we set F+ = 1 and Fx = 
two waveforms of a (1.4 + 10)Af(: 



0. In Fig. [T] we plot 
system to illustrate 



the effect of eccentricity on the waveform. In comparing 
the two waveforms, we find three major effects induced 
by eccentricity: (i) amplitude modulation; (ii) decreased 
duration of the signal; and (iii) increased signal ampli- 
tude. 





FIG. 1: Gravitational wave signals from a (1.4 -|- 10) Mq bi- 
nary generated using the x-model formalism. The waveforms 
start at /gw = 40 Hz and are terminated at the Schwarzschild 
ISCO frequency. The top panel shows a non-eccentric wave- 
form and in the bottom panel shows a signal with initial ec- 
centricity eo = 0.4. The sharp peaks are due to the increase 
in gravitational radiation that occurs during periastron pas- 
sage. The three major effects of eccentricity illustrated here 
are: (i) a decreases the duration of the waveform; (ii) wave- 
form amplitude modulation; and (iii) an overall increase in 
amplitude. 



III. DATA ANALYSIS 

LIGO's current searches for gravitational waves from 
neutron star and stellar-mass black hole binaries use 
the restricted stationary phase approximation to the 
Fourier transform of 3.5 FN circular waveforms (known 
as TaylorF2) [33|. These templates are given in the 
frequency domain by 

M/;M,ry) = A(M,,7)/-7/6e(/-/e)e»*(^^*^-''), (20) 

where A{M, rj) is the overall amplitude of the template 
at a canonical distance (typically 1 Mpc), 9 is the Heav- 
iside step function and fc is the upper cut-off frequency, 
given by the Schwarzschild innermost stable circular or- 
bit (ISCO) frequency, /isco = l/(6V67rAf). The am- 
plitude is included only at leading order, but the phase 
^'(/; M,ri) is computed to 3.5 FN order by 
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-V 




55 

—r] ) — 167rw'^ 
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V 756 

11583 231236 531 



/15293365 27145 
V 508 032 ^ 504 
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15 335 597 827 

3 048192 ^ ~l2 
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-77- 



4694 215 680 



127 825 
1296 



254 016 



1512 



756 



6 



21 



-V 



3085 



(7 + ln(4u)) 



(21) 



where v — {'wMff'^^ and 7 is Euler's constant [ssl. [39l|. 
To search for a signal in detector data, we construct the 
signal-to-noise ratio p from the template h and the cali- 
brated output of the detector s according to 

where the inner product is given by 



{s\h) = 2 df 



00 ~* 



S*{f)hif) + ~sif)h*if) 



Snif) 



f°° S*(f)h(f) 
s{f) denotes the Fourier transform of s{t), 



(23) 



(24) 



and S„{f) is the one-side noise power spectral density 
(PSD) of the detector. In the presence of Gaussian noise 
alone, the expectation value of is unity; large values of 
p indicate that a signal is present in the data. The tem- 
plate parameter tc gives the time of arrival of the si gnal , 
and can be searched over using a Fourier transform '40] . 
Similarly the coalescence phase (j)c can be searched over 
analytically [4l| : these parameters are termed extrinsic 
parameters. The masses of the signal are not known 
a priori^ however, and so a discrete bank of templates 
must be constructed and the detector data filtered 
against each template in the bank. Since the mass pa- 
rameters M, rj must be explicitly searched over, they are 
known as intrinsic parameters. The ambiguity function 
A{9) is defined with respect to the template parameters 
e^{tcAc,M,iq) as 



A{e) 



{s\h) 



^{s\s){h\h) 



(25) 



The maximization of Eq. ((25|) over the extrinsic template 
parameters defines the overlap 



Ois,h) 



tc.<t>c y/{s I s){h I h) 



— max {s \ h), 
tc,<f>c 



(26) 



|40l l42j , where the hat denotes normalization of the tem- 
plate amplitude. The fitting factor J-", quantifies loss in 
signal-to-noise ratio (SNR) due to the filtering of a signal 
with template that is not an optimal filter [ij] . The loss 
in event rate incurred by non-optimal filtering is propor- 
tional to 1 — J^^. For instance, a 3% decrease in gives 
a loss of approximately 10% in event rate. The value of 
is obtained by maximizing the overlap function over 
the intrinsic parameters of the template waveform 



= max max (s 

M, rj tc - 4>C 



h). 



(27) 



When constructing a bank, the templates are typically 
placed such that the loss in signal-to-noise ratio due to 
mismatch between a template with parameters inside the 
bank and the nearest bank grid point is no greater than 
3%; such a bank has a fitting factor of 0.97. 

In this paper, we use the hexagonal template place- 
ment algorithm from the LIGO Algorithm Library 
(LAL) [ij, to place a bank which runs over the mass 
range 1 < mi, TO2 < 34Mq subject to M < 35M0 at a 
fitting factor of 0.97. Template placement depends on the 
shape of the noise power spectrum of the detector [i^ . 
We model the current LIGO detectors by an analytic fit 
to the LIGO noise curve |4f| , given by 



5„(/) = 9.0 X 10 



-46 



(4.49j:) 



-56 



0.16a:"^-^^ + 0.32a;2 0.52 



Hz-\ (28) 



where x — //150IIz and / is frequency in Hz. Since the 
current LIGO detectors (Enhanced LIGO) have the same 
seismic isolation as the Initial LIGO detectors, they have 
essentially the same low frequency response. Therefore 
the Initial LIGO PSD is sufficient to model the Enhanced 
LIGO detectors for the purpose of this study. The Ad- 
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vLIGO PSD [43 is modeled by 
Snif) = 1-60 X 10-49 



300 

f 24 



f f 



-17 



V15Hz 

/300 a; 
I 90 



7 



f f 



-3.45 



V50Hz 
3.5 



110 



1.02 - 1.08 a;^ + 0.54 x*^ 
1 + 0.21x2 



Hz-i , (29) 



where x = //300Hz. The low- frequency cutoff /o for the 
templates is 40 Hz for LIGO and 10 Hz for AdvLIGO. 

To study the effectiveness of the template bank de- 
scribed above to capture signals from eccentric systems, 
we model eccentric binaries by the waveform given in 
Eq. (fT8|). which includes 3 PN conservative dynamics and 
2 PN radiative dynamics. The binary is chosen to be op- 
timally oriented [33|. We generate an eccentric signal 
and compute the overlap of the waveform against each 
template in the bank. The overlap is maximized over 
the t emp late bank to obtain the best match. Following 
Ref. [48j , we define the effective fitting factor T 



max \s 

h G bank 



h). 



(30) 



The function T provides a measure of how effective the 
template bank is in capturing a signal s that does not 
lie exactly in the template bank manifold. The distribu- 
tion of J- allows one to isolate the regions of the physi- 
cal parameter space where the mismatch is less than the 
minimal match of the bank. 



IV. RESULTS 

A. Initial and Enhanced LIGO 

In this section we present the results of our study for 
the current LIGO detectors. To understand the degree 
to which eccentricity affects the the x-model waveform, 
we investigated the overlap as a function of eccentricity 
using a circular (eo = 0) x-model waveform as the tem- 
plate and a x-model waveform with non-zero eccentricity 
as the signal. Fig. [5] shows the LIGO PSD overlap as 
a function of initial eccentricity for three different mass 
ratios. The overlaps are more sensitive to eccentricity 
for smaller mass systems, since the lower mass systems 
have longer durations which allow for the effects of ec- 
centricity to accumulate. We next investigate overlaps 
between a;-model signals and 2 PN TaylorT4 templates. 
In section III C[ we noted the similarities between the 
zero eccentricity limit of the x-model and the 2 PN Tay- 
lorT4 approximation. Comparisons are made with the 
TaylorT4 at 2 PN since that is the radiative PN order 
at which the z-model has been computed. In Fig. |31 
we plot the LIGO PSD overlap between 2PN TaylorT4 
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FIG. 2: The overlap (se = o |s) using the Initial LIGO PSD 
for three systems having total masses (1.4 -I- 1.4)Mq, (1.4 -|- 
5.O)M0 and (5.0 -I- 5.O)M0 . Here, Se = o denotes the circular 
x-model signal. The overlaps for the higher mass system are 
less affected by eccentricity because they have fewer cycles in 
the LIGO frequency band, and can thus be phase-matched 
against template waveforms via a relative shift in coalescence 
time without significant phase decoherence. 



waveforms and x-model waveforms of the same mass as 
a function of eccentricity. As the eccentricity approaches 
zero, the overlaps between the T4 and the x-model ap- 
proach unity for all three signal masses, consistent with 
the zero eccentricity limit of the a;-model in Sec. Ill CI 

Fig. U shows overlaps for Initial LIGO between 2 PN x- 
model waveforms and 3.5 PN TaylorF2 templates. The 
overlap is worse for these waveforms, as would be ex- 
pected. At Co = the reduction in overlap is due to the 
difference between the PN models used to generate wave- 
forms used as signals and templates. As eo increases, 
the figure shows the additional loss in overlap due to 
the effect of eccentricity. The fact that the overlaps are 
not unity at cq = is due to the differences between 
the PN formulation of the time-domain x-model signal, 
which contains radiative effects to 2 PN order, and the 
frequency-domain TaylorF2 template, which contains ra- 
diative effects to 3.5 PN order. Although higher-order 
post-Newtonian effects are more significant (at a given 
frequency) for higher-mass binaries, more massive bina- 
ries coalesce at lower frequencies. At cq = 0, the lowest 
overlap is observed in Af = 2.8 Mq signal. Although 
higher order post-Newtonian corrections are small for 
this system, they have many waveform cycles over which 
to accumulate, resulting in a reduction of overlap. Note, 
however, that in calculating the overlaps shown in Figs. 
|3]and[5]the intrinsic parameters are kept fixed. Below we 
show the effect of maximizing the overlap over a bank, 
as is the case in a real search. As eg increases, eccentric- 
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FIG. 3: Overlaps using the LIGO PSD between the 2 PN 
x-model waveform and the 2 PN Taylor T4 waveform of the 
same total mass for masses (1.4+1.4) Mq, (1.4 + 5.0) Mq and 
(5.0 + 5.0) Mq. The overlaps approach one as the eccentricity 
goes to zero, which indicates that the two models are in good 
agree in the limit eo 0. 



ity causes additional higher-order post-Newtonian cor- 
rections to the signal waveform, but also causes the du- 
ration of the waveform to decrease. In the M = 6.8 Mq 
and M = 10 Mq waveforms, competition between these 
effects causes the overlap curves to cross at ep — 0.05, 
with slightly larger overlaps in the higher-mass signal at 
larger eccentricity. 



M = 2.8 A/q 

i\/ = 6.4 M@ 

- ■ - ■ - A/ = 10 A/r 




FIG. 4; Initial LIGO overlaps between the 2 PN a::-model and 
the 3.5 PN TaylorF2 waveforms as a function of initial eccen- 
tricity for (1.4 + 1.4) Mq, (1.4 + 5.0) Mq and (5.0 + 5.0) Mq 
systems. Notice that the eo = overlaps are not unity. This 
is due to differences in the PN order of the waveforms 



We now explore the effectiveness of the template bank 
described in Sec.|nT]at detecting eccentric signals. To ex- 
plore this for a range of signal masses and eccentricities, 
we generate a rectangular grid of 10 026 2:-model signals 
with masses in the range IMq < mi, m2 < UMq sub- 
ject to 2 < M /Mq < 15 and initial eccentricities in the 
range < eo < 0.4. Each signal in this grid is filtered 
through the 14 863 templates in the 3.5 PN TaylorF2 Ini- 
tial LIGO template bank. The effective fitting factor 
of the bank is computed for each signal. The signal space 
is three dimensional as it depends on mi, m2, and eo, and 
hence ^ is a function of three variables. To visualize the 
results of the bank simulation, we plot as a function 
of the total mass M and eccentricity eg of the injected 
x-model signal. In our signal grid, several component 
masses have the same total mass and so in Fig. [SJ we 
plot the highest, lowest, and mean T for each value of 
M. Notice that by maximizing over the template banks, 
we have gained significantly in overlap for the lower mass 
systems (cf. Fig. 2]). In the worst case, the 3.5 PN Tay- 
lorF2 bank can achieve effective fitting factors > 0.96 for 
all x-model signals which have a total mass M > 3Mq 
and an eccentricity cq < 0.05 at 40 Hz. The signals 
with the worst values of J- are those which are matched 
against templates which lie near the edges of the tem- 
plate bank. To separate the effect of eccentricity from 
effects related to the difference in PN order between the 
signal and template we calculate the relative difference 
between the effective fitting factor for signals with zero 
eccentricity and those with non-zero eccentricity, given 
by 



A^(eo) 



-^(eo = 0) 



(31) 



Fig. [7] shows for the highest, lowest, and average 
values of T. These distributions illustrate that the loss 
in overlap results predominantly from eccentricity at 
higher masses, but at lower masses is the accumulation of 
higher-order PN effects (included in the templates, but 
not in the signals) over the longer waveforms. We there- 
fore conclude that the current LIGO template bank is 
effective at capturing signals with small residual eccen- 
tricities eo < 0.05 at 40 Hz for all masses 2 < M < 15M0_. 
For eccentricities higher than eo ~ 0.05, the loss in T 
becomes more significant, although for the higher mass 
systems this loss does not occur until eo > 0.1. 



B. Results for Advanced LIGO 

Finally, we investigate how eccentricity effects the de- 
tection efficiency of TaylorF2 templates for AdvLIGO. 
Although it is unlikely that TaylorF2 templates will be 
used for AdvLIGO searches, this study is illustrative of 
the effect of eccentricity on circular templates with the 
AdvLIGO noise curve. The lower frequency cut-off for 
the AdvLIGO PSD requires that we start our waveforms 
at 10 Hz instead of 40 Hz. Fig. [S] shows the overlap 
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FIG. 5: Overlaps between the 2 PN a;-model eccentric wave- 
forms and the 3.5 PN TaylorF2 waveforms using the Ad- 
vLIGO PSD. Reduced overlaps result from the increase in the 
amount of time spent in the detector's band. The increase in 
the duration leads to larger phase de-coherence and intensifies 
the effects due to differences in PN order. Additionally, the 
difference in the PN order of the signals and the templates 
reduces the overlap even at eo = 0. 



between the 3.5 PN TaylorF2 waveforms and the eccen- 
tric x-model waveforms for the AdvLIGO noise curve. 
Since the waveforms have more detectable cycles in Ad- 
vLIGO, a larger phase difference accumulates leading to 
lower overlaps than for Initial LIGO. The effect of the 
post-Newtonian differences between the waveform and 
templates is also apparent in Fig. [5l as in Fig. S) 

Our analysis of the AdvLIGO template bank results 
proceeds in a similar way to that for the current de- 
tectors. Due to the improved low-frequency sensitivity 
of AdvLIGO, the size of the template bank increases 
to 136 000, greatly increasing the computational cost of 
the template bank simulation. As a result, we decrease 
the resolution of the signal grid to inject only 280 sig- 
nals. This is sufficient to probe the effect of eccentricity 
on AdvLIGO searches, however. Fig. [5] shows the best, 
worst and mean values of f . Notice that J- is very poor 
for the long, low-mass signals, and the high-mass sig- 
nals (where higher order PN effects dominate) even at 
zero eccentricity. This is due to the difference in PN or- 
der between the signals and the templates. In Fig. [3] we 
plot for the highest, lowest, and average values of F 
for AdvLIGO. In the AdvLIGO case, the effect of small 
residual eccentricity can be significant for low mass sys- 
tems {M < AMq). Further studies using both higher 



PN-order eccentric waveforms and the actual AdvLIGO 
search templates are needed to confirm these findings. 
V. CONCLUSION 

In this paper we studied the effect of eccentricity on 
detection searches for gravitational waves from compact 
binary coalescence using the current and Advanced LIGO 
detectors. To construct eccentric binary waveforms, we 
used the NR-calibrated 2 PN x-model formalism pro- 
posed by Hinder et al. [2^. We considered the zero ec- 
centricity limit of the x-model and compared it with the 
TaylorT4 approximant; analytic and overlap calculations 
confirm that the phase evolutions of these two models 
agree in the eo = limit. Using x-model waveforms 
as our signal family, we performed template bank sim- 
ulations using the 3.5 PN TaylorF2 waveforms and the 
hexagonal placement algorithm implemented in L AL p5| . 
The 3.5 PN TaylorF2 template bank was found to give 
fitting factors of 0.96 or better for systems with total 
mass M > 3Mq and eccentricity at 40 Hz of eo < 0.05. 
By separating the loss in fitting factor due to differences 
in higher-order circular PN corrections and the effect 
of eccentricity, we conclude that current LIGO searches 
are sensitive to binaries with small residual eccentricities 
when the waves enter the sensitive band of the detec- 
tor. However for eccentricities eo 0.1, significant losses 
in sensitivity will be observed. For AdvLIGO, our re- 
sults were dominated by differences in the PN order of 
the signals and templates and so more investigation with 
higher-order PN eccentric signals is needed. However, 
our study suggests that AdvLIGO is even more sensitive 
to the effect of residual eccentricity in compact binary 
inspiral, as would be expected. 

In this study, we did not consider the ability of the 
LIGO detectors to measure eccentricity, however our 
results suggest that for eo > 0.1, LIGO is sensitive 
to eccentricity, even in the current detectors. Careful 
study of parameter estimation and implementation of a 
search for eccentric systems with Enhanced LIGO or Ad- 
vLIGO would allow the formation mechanisms proposed 
in Refs. [H, [l3| to be explored. 
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